On the kinetic energy of unitary Fermi gas in a harmonic trap 
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We have considered the orbital-free approximation of the kinetic energy func- 
tional to investigate the zero temperature properties of dilute harmonically trapped 
two component Fermi gas at unitarity. It is shown that our approach provides 
a realible and inexpensive method to study superfluid strongly interacting dilute 
Fermi gases. 
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The modern density functional theory (DFT) is based on the Kohn-Sham (KS) ap- 
proach [1], where the noninteracting kinetic energy, T, is calculated in terms of the KS 
orbitals, although Ref.[2] proved the basic existence of functional T(p), where p is the density, 
/ p{r)d^r = N, and is the particle number. Since the numerical cost of self-consistently 
determining N orbitals rapidly increases for large N, the accurate density-functional approx- 
imation to the kinetic energy in terms of the density would reduce dramatically complex- 
ity of the DFT calculations (here we note the superfiuid extension of the DFT given in 
Refs.[3-5]). There are various noninteracting kinetic energy functionals proposed in litera- 
ture, for chemical applications see, for example, Refs.[6-21] (for nonlocal density functionals 
see Refs. [22-26]). For applications of the DFT to the nuclear structure physics see Ref. [27], 
web site, constructed for the universal nuclear energy density functional (UNEDF) collabo- 



ration, http://unedf.org| and references therein. 



The kinetic energy functional can be written as 
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T[p] = ^ J r{p{f))d\ (1) 
where the Kirgnitz semi classical expansion for the kinetic energy density [28-31] 

r(p(f)) = TTF{p) + lrw{p) + ..., (2) 

where 

TTFip) = liSnY'p'^', (3) 

is the well-known Thomas- Fermi (TF) kinetic energy density and Tw[p) is the Weizsacker 
kinetic energy density [12] 

rwip) = [^p'/\r-)r. (4) 

The semiclassical expansion (2) has to be considered as an asymptotic expansion [11]. While 
ttf is exact for the uniform gas model, Tw{p) is considered exact in the limit of rapidly 
varying density p [13,14]. 

We present in this paper the simple linear combinations 

T = Xi{N)ttf + HN)tw (5) 

with Aj determined empirically from the asymptotic region and for getting good energy of 
few fermion systems at unitarity. A N-dependence of A2 corresponds to the effective mass 
N-dependence [32]. 



There has been a lot of interest in systems of fermions at the unitarity [32-35] (when 
the scattering length diverges, the Bertsch many-body problem [36]). While Refs. [37-40] 
consider homogeneous systems, Refs. [41-45] present ab initio calculations of the properties 
of trapped fermionic atoms. 

In Refs. [46,47] and later in Refs. [48-63] the dynamics of strongly interacting trapped di- 
lute Fermi gases (dilute in the sense that the range of interatomic potential is small compared 
with inter-particle spacing) consisting of a 50-50 mixture of two different states is investi- 
gated in the single equation approach to the time-dependent density-functional theory, using 
Ai = 1 and a constant A2 approximations. 

For the stationary case Eq.(5) leads to the following DFT equation 

-X2{N)—V'^ + Kxt* + V^c^ = (6) 

where Vxc{r) = [^^^]p=p{r)) e(p) is the ground-state energy per particle of the homogeneous 
system, p{f) =\ ^(r) p and /i is the chemical potential. For the remainder of this paper we 
will consider fermion systems at unitarity in a spherical harmonic trap 

Vext{r) = . (7) 

The ground state energy is given by the minimum of the energy functional 

J[m] = A2(iV)|- J Tw{p)dh + j Ve^tpdh + j e{p)pd\ (8) 

where e(p) = (Ai(A^) + P)3h'^kp/{10m), kp ^ (37r^p)^/^, p = and the universal param- 
eter P [33] is estimated tohe P — —0.56 [39]. 

To test the accuracy of approximation (5), we write the radial nonlinear equation (6) in 
asymptotic region 

The regular as r — >• 00 solution of Eq.(9) can be written as ^^Wi^^i/ 2+1 /ii^)-, where W is the 
regular as r — >• 00 Whittaker function, z = {mui/{h^J \2{N)))r'^ and n = p,/ {2?iu!^ X2{N)), 
therefore 

p-e-V^-^/^ (10) 

Since, in the limit of large r the Hartree-Fock density, PijF(r) is proportional to the square 
of the last occupied state 

pHF{r) - e-"^'*^'■'/V'^/('''^)-^ (11) 



we expect that limjv-*oo '^2(-^) — 1- We also expect that \imN_^^ Xi{N) — > 1, since the TF 
kinetic energy density becomes exact in the large N limit. Note, that in the large r region 
the density is not slowly varying. For the kinetic energy functional in one spatial dimension 
Ref.[64] derived an upper bound 

< r{p) ><< Tw{p) > + < ttf{p) > ■ (12) 

Although, the question whether the upper bound, Eq.(12), holds in the three dimensional 
case is still open question [65] , our numerical results are in strong support of the inequality 
(12), see Fig. 1. 

To study the effectiveness of the approximations (5) we calculate lower bounds to the 
ground state energy. Eg. (8), using the results of Rcf.[62] 



B(-. = 5s.^^/A.w + 5i^3MW. (13) 

To calculate upper bounds, , we employing Fetter's trial functions [66] 

pV2(^-3^i,(l_(l_^)(rf^)2)i/M^ (14) 

where d and q are the variational parameters and b is the normalization constant, to minimize 
the functional J, Eq. (8). 

Following Ref.[24], we propose the following approximation for the kinetic energy density, 
X2 — 1 and 

X,(N) = {l-j^){l-^), (15, 

where c and 7 are fixed by a least squares fit to the fixed-node diffusion Monte Carlo data 
[42], c = 1.46832, 7 = 0.78383. It is clear from Table I that the fit is a very accurate. As for 
the lower and the upper bounds, they provide the actual solution of equation (6), 
E = (£;(+) + £;(-))/2 within ±S accuracy, with S < 1%. 

Recently, Refs. [58-60] have considered nonhnear equation, which for stationary case cor- 
responds to the following approximation of the kinetic energy density 

1 , , 

T = -Tw + TtF- (16) 

Figures 2 and 3 show comparison between energy calculations using approximations (15), 
(16), (12) and the Thomas- Fermi approximation. It indicates that (i) there is a very good 
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agreement between calculations using (15) and (16) for > 10, (ii) the difference between 
all four approximations is negligible for N/2 > 10^. 

The T= {1/4:)tw + ttf approximation has incorrect asymptotic behavior for large r and 

N. The success of this approximation (see figures 2 and 3), however, indicates that the 
large r behavior have little or no consequence in the calculation of ground state energy of N 
fermion systems at unitarity in a spherical harmonic trap. 
In conclusion, we summarize the main points of this paper. 

(i) We have considered the orbital free approximation of the kinetic energy functional, 
proposed for the first time by Achariya et al. [19], to investigate N-fermion systems at 
unitarity consisting of 50-50 mixture of two different states and confined in a spherical 
harmonic trap. 

(ii) We found that our analytical lower bound, Eq.(13), describes the ground state energy 
with a very good accuracy, providing an easy and simple quantitative tool for trapped Fermi 
gases, without relying on complex and extensive computations. 

A.L.Z thanks N.J. Giordano and W.L. Fornes for providing the opportunity to finish this 
work. 
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Table I. The energies E^-\ Eq. (), E^+\E = {E^'^ + E^+^)/2 and the energy calculated 
within the fixed-node diffusion Monte Carlo method, E^'^ [42] , all in units of ku; for < 30 
[see the text for further details] . 
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Fig.l. Ground state energy per particle of few-fermion systems at unitarity in a spherical 
harmonic trap in units of hu as a function of number of atoms N . The solid line represents 
T — ttf -\- Tw approximation. The circular dots indicate results of Ref.[42]. 
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Fig. 2. Ground state energy per particle in units of huj as a function of number of atoms 
N . The solid line, the circular dots, the dashed line and the dashed-dotted line represent 
results calculated using t = tw + ^i{N)ttf, t = (1/4)tvk + ttf, t = tw + 'Ttf and r = ttf 
approximation, respectively. 
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Fig. 3. Same as in Fig. 2. 
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